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Abstract 

I We illustrate the use of machine learning techniques to analyze, structure, maintain, 

q ■ and evolve a large online corpus of academic literature. An emerging field of research 

can be identified as part of an existing corpus, permitting the implementation of a 
Q ' more coherent community structure for its practitioners. 
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1 Background 



The arXiv 1 is an automated repository of over 250,000 full-text research articles 2 in physics 
and related disciplines, going back over a decade and growing at a rate of 40,000 new sub- 
missions per year. It serves over 10 million requests per month [Ginsparg, 2003], including 
tens of thousands of search queries per day, and over 20 million full text downloads during 
calendar year '02. It is a significant example of a Web-based service that has changed the 
practice of research in a major scientific discipline. It provides nearly comprehensive cover- 
age of large areas of physics, and serves as an on-line seminar system for those areas. It also 
provides a significant resource for model building and algorithmic experiments in mapping 
scholarly domains. Together with the SLAC SPIRES-HEP database 3 , it provides a public 
resource of full-text articles and associated citation tree of many millions of links, with a 
q ! focused disciplinary coverage, and rich usage data. 

In what follows, we will use arXiv data to illustrate how machine learning methods can be 
used to analyze, structure, maintain, and evolve a large online corpus of academic literature. 
The specific application will be to train a support vector machine text classifier to extract 
an emerging research area from a larger-scale resource. The automated detection of such 
subunits can play an important role in disentangling other sub-networks and associated sub- 
communities from the global network. Our notion of "mapping" here is in the mathematical 
sense of associating attributes to objects, rather than in the sense of visualization tools. 
While the former underlies the latter, we expect the two will increasingly go hand-in-hand 
(for a recent review, see [Bonier et al., 2003]). 

*Presenter of invited talk at Arthur M. Sackler Colloquium on "Mapping Knowledge Domains", 9-11 
May 2003, Beckman Center of the National Academy of Sciences, Irvine, CA 
1 See http://arXiv.org/ . For general background, see [Ginsparg, 2001]. 
2 as of mid-Oct 2003 

3 The Stanford Linear Accelerator Center SPIRES-HEP database has comprehensively catalogued the 
High Energy Particle Physics (HEP) literature online since 1974, and indexes more than 500,000 high-energy 
physics related articles including their full citation tree (see [O'Connell, 2002]). 
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Figure 1: Representing text as a feature vector. 



2 Text Classification 

The goal of text classification is the automatic assignment of documents to a fixed number 
of semantic categories. In the "multi-label" setting, each document can be in zero or one or 
more categories. Efficient automated techniques are essential to avoid tedious and expensive 
manual category assigment for large document sets. A "knowledge engineering" approach, 
involving hand-crafting accurate text classification rules, is surprisingly difficult and time- 
consuming [Hayes and Weinstein, 1990]. We therefore take a machine learning approach to 
generating text classification rules automatically from examples. 

The machine learning approach can be phrased as a supervised learning problem. The 
learning task is represented by the training sample S n 

(x 1 ,y 1 ),(x 2 ,y2),---,(xn,y n ) (1) 

of size n documents, where Xi represents the document content. In the multi-label setting, 
each category label is treated as a separate binary classification problem. For each such 
binary task, G { — 1, +1} indicates whether a document belongs to a particular class. The 
task of the learning algorithm C is to find a decision rule he '■ x — > { — 1, +1} based on S n 
that classifies new documents x as accurately as possible. 

Documents need to be transformed into a representation suitable for the learning algo- 
rithm and the classification task. Information Retrieval research suggests that words work 
well as representation units, and that for many tasks their ordering can be ignored without 
losing too much information. This type of representation is commonly called the "bag-of- 
words" model, an attribute-value representation of text. Each text document is represented 
by a vector in the lexicon space, i.e., by a "term frequency" feature vector TF(iUi,x), with 
component values equal to the number of times each distinct word Wi in the corpus occurs 
in the document x. Fig. 1 shows an example feature vector for a particular document. 

This basic representation is ordinarily refined in a few ways: 

TFxIDF Weighting: Scaling the components of the feature vector with their inverse doc- 
ument frequency IDF(wj) [Salton and Buckley, 1988] often leads to improved perfor- 
mance. In general, IDF(wj) is some decreasing function of the word frequency DF(iUj), 
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equal to the number of documents in the corpus which contain the word Wi. For 
example, 



IDF(^) = log 



( 



DF(iUi) 



n 



) 



(2) 



where n is the total number of documents. Intuitively, the inverse document frequency 
assumes that rarer terms have more significance for classification purposes, and hence 
gives them greater weight. To compensate for the effect of different document lengths, 
each document feature vector Xi is normalized to unit length: = 1. 

Stemming: Instead of treating each occurrence form of a word as a different feature, stem- 
ming is used to project the different forms of a word onto a single feature, the word 
stem, by removing inflection information [Porter, 1980]. For example "computes", 
"computing" , and "computer" are all mapped to the same stem "comput" . The terms 
"word" and "word stem" will be used synonymously in the following. 

Stopword Removal: For many classification tasks, common words like "the", "and", or 
"he" do not help discriminate between document classes. Stopword removal describes 
the process of eliminating such words from the document by matching against a pre- 
defined list of stop-words. We use a standard stoplist of roughly 300 words. 



SVMs [Vapnik, 1998] were developed by V. Vapnik et al. based on the structural risk mini- 
mization principle from statistical learning theory. They have proven to be a highly effective 
method for learning text classification rules, achieving state-of-the-art performance on a 
broad range of tasks [Joachims, 1998, Dumais et al., 1998]. Two main advantages of using 
SVMs for text classification lie in their ability to handle the high dimensional feature spaces 
arising from the bag-of-words representation. From a statistical perspective, they are robust 
to overfitting and are well suited for the statistical properties of text. From a computational 
perspective, they can be trained efficiently despite the large number of features. A detailed 
overview of the SVM approach to text classification, with more details on the notation used 
below, is given in [Joachims, 2002]. 

In their basic form, SVMs learn linear decision rules 



described by a weight vector w and a threshold b, from an input sample of n training examples 
S n = Ui)) ' ' ' ) {%ni Un))) %i £ IR^, Ui € {— 1, +1}. For a linearly separable S n , the SVM 
finds the hyperplane with maximum Euclidean distance to the closest training examples. 
This distance is called the margin 8, as depicted in Fig. 2. Geometrically, the hyperplane 
is defined by its normal vector, w, and its distance from the origin, —b. For non-separable 
training sets, the amount of training error is measured using slack variables £j. 

Computing the position of the hyperplane is equivalent to solving the following convex 
quadratic optimization problem [Vapnik, 1998]: 



3 Support Vector Machines 



h(x) = sgn{w ■ x + b} , 



(3) 
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Figure 2: A binary classification problem (+ vs. — ) in two dimensions. The hyperplane h* 
separates positive and negative training examples with maximum margin 5. The examples 
closest to the hyperplane are called support vectors (marked with circles). 



Optimization Problem 1 (SVM (primal)) 

minimize: V(w, b, £) — -w ■ w + C (4) 

2 i=i 

subj. to: ^ =1 :y i [w-x i + b\>l-li i (5) 

K=i :^>0 (6) 

The margin of the resulting hyperplane is 5 — 

The constraints (5) require that all training examples are classified correctly up to some 
slack £j. If a training example lies on the "wrong" side of the hyperplane, we have the 
corresponding £j > 1, and thus X^Li C« i s an upper bound on the number of training errors. 
The factor C in (4) is a parameter that allows trading off training error vs. model complexity. 
The optimal value of this parameter depends on the particular classification task and must be 
chosen via cross-validation or by some other model selection strategy. For text classification, 
however, the default value of C — l/maxj||xj|| 2 = 1 has proven to be effective over a large 
range of tasks [Joachims, 2002]. 

OP1 has an equivalent dual formulation: 

Optimization Problem 2 (SVM (dual)) 

n -y n n 

maximize: W(d) =^ a, - - J2 ViVj^A^ ' fy) ( 7 ) 
i=i z i=ij=\ 



subj. to: Vi a i = ( 8 ) 
i=i 

Vi G [l..n] : < Qi < C (9) 
From the solution of the dual, the classification rule solution can be constructed as 

n 

w = ^2 aiyA and b = y usv - w-x usv , (10) 
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where (x usv , y usv ) is some training example with < a usv < C. For the experiments in this 
paper, SVM Llght [Joachims, 2002] is used for solving the dual optimization problem 4 . More 
detailed introductions to SVMs can be found in [Burges, 1998, Cristianini and Shawe- Taylor, 2000]. 

An alternative to the approach here would be to use Latent Semantic Analysis (LSA/LSI) 
[Landauer and Dumais, 1997] to generate the feature vectors. LSA can potentially give bet- 
ter recall by capturing partial synonymy in the word representation, i.e., by bundling related 
words into a single feature. The reduction in dimension can also be computationally efficient, 
facilitating capture of the full text content of papers, including their citation information. 
On the other hand, the SVM already scales well to a large feature space, and performs well 
at determining relevant content through suitable choices of weights. On a sufficiently large 
training set, the SVM might even benefit from access to fine distinctions between words 
potentially obscured by the LSA bundling. It is thus an open question, well worth further 
investigation, as to whether LSA applied to full text could improve performance of the SVM 
in this application, without loss of computational efficiency Generative models for text 
classification provide yet another alternative approach, in which each document is viewed as 
a mixture of topics, as in the statistical inference algorithm for Latent Dirichlet Allocation 
used in [Griffiths and Steyvers, 2003]. This approach can in principle provide significant ad- 
ditional information about document content, but does not scale well to a large document 
corpus, so the real-time applications intended here would not yet be computationally feasible 
in that framework. 

4 arXiv Benchmarks 

Before using the machine learning framework to identify new subject area content, we first 
assessed its performance on the existing (author-provided) category classifications. Roughly 
180,000 titles and abstracts were fed to model building software which constructed a lexi- 
con of roughly 100,000 distinct words and produced training files containing the TDxIDF 
document vectors for SVM Llght . (While the SVM machinery could easily be used to analyze 
the full document content, previous experiments [Joachims, 2002] suggest that well- written 
titles and abstracts provide a highly focused characterization of content at least as effective 
for our document classification purposes.) The set of support vectors and weight parameters 
output by SVM Llght was converted into a form specific to the linear SVM, eq. (3): a weight 
vector w c and a threshold b c , where c is an index over the categories. 

As seen in Fig. 3, the success of the SVM in classifying documents improves as the 
size of a category increases. The SVM is remarkably successful at identifying documents in 
large (> 10, 000 documents) categories and less successful on smaller subject areas (< 500 
documents). A cutoff was imposed to exclude subject areas with fewer than 100 documents. 5 

In experiments with small categories (N < 1000), the SVM consistently chose thresholds 
that resulted in unacceptably low recall (i.e., missed documents relevant to the category). 
This is in part because the null hypothesis, that no documents are in the category, describes 

4 SVM Llght is available at http://svmlight.joachims.org/ 

5 Some of the smaller subject areas are known to be less topically focused, so the difficulty in recall, based 
solely on title/abstract terminology, was expected. 
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Figure 3: Recall and Precision as functions of category size for 78 arXiv major categories 
and minor subject classes. 2/3 of the sample was used as a training set and 1/3 as a test 
set. The four largest categories, each with over 30,000 documents are cond-mat, astro-ph, 
hep-th, and hep-ph. 



the data well for a small category: e.g., if only 1000 of 100,000 documents are in a category, 
the null hypothesis has a 99% accuracy. 6 To counteract the bias against small categories, 
10% of the training data was used as a validation set to fit a sigmoid probability distribution 
l/(l+exp(Af+B)) [Piatt, 1999]. This converts the dot product Xi ■ w c between document 
vector and weight vector into a probability P(i e c | Xi) that document i is a member of 
category c, given its feature vector X{. Use of the probability in place of the uncalibrated 
signed distance from the hyperplane output by the SVM permitted greater levels of recall 
for small categories. 

Other experiments showed that the use of TFxIDF weighting as in eq. (2) improved 
accuracy consistently over pure TF weighting, so TFxIDF weighting was used in the exper- 
iments to follow. We also used a document frequency threshhold to exclude rare words from 
the lexicon, but found little difference in accuracy between using a document occurrence 
threshold of two and five. 7 Increasing the weight of title words with respect to abstract 
words, on the other hand, consistently worsened accuracy, indicating that words in a well- 
written abstract contain as much or more content classification import as words in the title. 
Changes in the word tokenization and stemming algorithms did not have a significant impact 
on overall accuracy. The default value of C — 1 in eq. (9) was preferred. 

6 "Accuracy" is the percentage of documents correctly classified. We also employ other common termi- 
nology from information retrieval: "precision" is the fraction of those documents retrieved that relevant to 
a query, and "recall" is the fraction of all relevant documents in the collection that are retrieved. 

7 Words that appeared in fewer than two documents constituted roughly 50% of the lexicon, and those 
that appeared in fewer than five documents roughly 70%. Ignoring rare and consequently uninformative 
words hence reduces the computational needs. 
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5 q-bio Extraction 



There has been recent anecdotal evidence of an intellectual trend among physicists towards 
work in biology, ranging from biomolecules, molecular pathways and networks, gene expres- 
sion, cellular and multicellular systems to population dynamics and evolution. This work 
has appeared in separate parts of the archive, particularly under "Soft Condensed Mat- 
ter" , "Statistical Mechanics" , "Disordered Systems and Neural Networks" , "Biophysics, and 
"Adaptive and Self- Organizing Systems" (abbreviated cond-mat.soft, cond-mat.stat-mech, 
cond-mat.dis-nn, physics. bio-ph, and nlin.AO). A more coherent forum for the exchange 
of these ideas was requested, under the nomenclature "Quantitative Biology" (abbreviated 
"q-bio"). 

To identify first whether there was indeed a real trend to nurture and amplify, and to 
create a training set, volunteers were enlisted to identify the q-bio content from the above 
subject areas in which it was most highly focused. Of 5565 such articles received from Jan 
2002 through Mar 2003, 466 (8.4%) were found to have one of the above biological topics as 
its primary focus. The total number of distinct words in these titles, abstracts, plus author 
names, was 23558, of which 7984 were above the DF = 2 document frequency threshold. 
(Author names were included in the analysis since they have potential "semantic" content 
in this context, i.e., are potentially useful for document classification. The SVM algorithm 
will automatically determine whether or not to use the information by choosing suitable 
weights.) 

A data-cleaning procedure was employed, in which SVM Llght was first run with C = 
10. We recall from eqs. (4) and (9) that larger C penalizes training errors and requires 
larger ct's to fit the data. Inspecting the "outlier" documents [Guyon et al., 1996] with the 
largest \ai\ then permitted manual cleaning of the training set. 10 were moved into q-bio, 
and 15 moved out, for a net movement to 461 q-bio (8.3%) of the 5565 total. Some of 
the others flagged involved word confusions, e.g., "genetic algorithms" typically involved 
programming rather than biological applications. Other "q-bio" words with frequent non- 
biological senses were "epidemic" (used for rumor propagation), "evolution" (used also for 
dynamics of sandpiles), "survival probabilities", and extinction. "Scale-free networks" were 
sometimes used for biological applications, and sometimes not. To help resolve some of these 
ambiguities, the vocabulary was enlarged to include a list of most frequently used two-word 
phrases with semantic content different from their constituent words . 

With a training set fully representative of the categories in question, it was then possible 
to run the classifier on the entirety of the same subject area content received from 1992 
through Apr 2003, a total of 28,830 documents. The results are shown in Fig. 4. A total 
of 1649 q-bio documents was identified, and the trend towards an increasing percentage of 
q-bio activity among these arXiv users is evident: individual authors can be tracked as they 
migrate into the domain. Visibility of the new domain can be further enhanced by referring 
submitters in real time, via the automated classifier running behind the submission interface. 

Some components of the weight vector generated by the SVM for this training set are 
shown in Fig. 5. Since the distance of a document to the classifying hyperplane is determined 
by taking the dot product with the normal vector, its component values can be interpreted 
as the classifying weight for the associated words. The approach here illustrates one of the 
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Figure 4: The number of submissions per year from 1992 through April 2003 in particular 
subsets of arXiv subject areas of cond-mat, physics, and nlin most likely to have "quantitative 
biology" content. The percentage of q-bio content in these areas grew from roughly 1% to 
nearly 10% during this timeframe, suggesting a change in intellectual activity among arXiv- 
using members of these communities. 
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protein +8.57 
dna +7.08 
biological +5.06 
neuron +5.05 
rna +3.30 
gene +3.21 
mutation +3.15 
population +3.11 
epidemic +3.05 
biology +3.02 
disease +2.93 

cell +2.90 
neural +2.89 
brain +2.83 
ecosystem +2.56 

tissue +2.52 
sequence +2.51 
genetic +2.51 
bacterial +2.48 
blood +2.43 
genome +2.37 
peptide +2.37 
infection +2.34 



forward +0.00 
minimalist +0.00 

region +0.00 
confinement +0.00 
implies +0.00 
96 +0.00 
y_togashi +0.00 
n_wingreen +0.00 
mean_free +0.00 
narrower +0.00 
shot -0.00 
repton -0.00 
kyoto -0.00 
regular -0.00 
generalisation -0.00 
d^saakian -0.00 
conformity -0.00 

aware -0.00 
evemthough -0.00 
practitioner -0.00 
permittivity -0.00 



point -1.04 
equation -1.05 
boundary -1.06 
social -1.06 
n -1.09 
relaxation -1.14 
fluid -1.15 
indian -1.15 
spin -1.17 
spin_glass -1.17 
traffic -1.18 
system -1.30 
polymer -1.33 

class -1.35 
emerge -1.36 
gradient -1.39 
quantum -1.43 
surface -1.43 
synchronization -1.45 
market -1.47 
particle -1.52 
polyelectrolyte -1.53 
world -1.57 



Figure 5: Shown above are the most positive, a few intermediate, and most negative com- 
ponents of the q-bio classifying weight vector. 
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major lessons of the past decade: the surprising power of simple algorithms operating on 
large datasets. 

While implemented as a passive dissemination system, the arXiv has also played a social 
engineering role, with active research users developing an affinity to the system and adjusting 
their behavior accordingly. They scan new submissions on a daily basis, assume others in 
their field do so and are consequently aware of anything relevant that has appeared there 
(while anything that doesn't may as well not exist), and use it to stake intellectual priority 
claims in advance of journal publication. We see further that machine learning tools can 
characterize a subdomain and thereby help accelerate its growth, via the interaction of an 
information resource with the social system of its practitioners. 8 

Postscript: The q-bio extraction described in section 5 was not just a thought experiment, 
but a prelude to an engineering experiment. The new category went on-line in mid-September 
2003, 9 at which time past submitters flagged by the system were asked to confirm the q-bio 
content of their submissions, and to send future submissions to the new category. The activity 
levels at the outset corresponded precisely to the predictions of the SVM text classifier, and 
later began to show indications of growth catalyzed by the public existence of the new 
category. 
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